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Abstract 

Two problems in the search of metric characteristics on weighted 
' undirected graphs with non-negative edge weights are being considered. 

The first problem: a weighted undirected graph with non-negative edge 
weight is given. The radius, diameter and at least one center and one 
pair of peripheral vertices of the graph are to be found. In the second 
Q-j problem we have additionally calculated the distances matrix. For the 

^ . problems being considered, we proposed fast search algorithms which 

use only small fraction of graph's vertices for the search of the metric 
characteristics. The proposed algorithms have been compared to other 
popular methods of solving problems considered on various inputs. 
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Q ! Introduction 

Zj . The metric (or numerical) characteristics of a graph are parameters 

defined by the shortest paths between vertices: center, radius and diameter. 
The diameter of a graph is the longest of the shortest paths between all pairs 
of vertices in the graph. The center is a vertex with the smallest maximum 
\ distance to any other vertex in the graph. This smallest distance is the 

t^- 1 radius of the graph. 

These characteristics are essential properties of graphs, they represent 
ON I valuable information about structures corresponding to the graphs. For 

1 example, in the facility location problem, the diameter can represent the 

maximum length of a waiting period, center and radius " — are the best 
location of the facility and the maximum time needed to reach this facility. 
In cases where graphs correspond to computer networks, the diameter can 
represent a connection speed between the two slowest nodes of the network 
and the radius can display the value of a delay between a server and the 
slowest node etc. 

In general, the method of searching for these characteristics is to find 
the shortest paths between all pairs of vertices of a given graph and to 
determine distances and vertices which satisfy the definitions of the center, 
radius and diameter. Known algorithms performing this task on weighted 
graphs with n vertices have a complexity from 0(Cn 2 ' 376 ) for graphs with 
integer edge weights less than C [11] and up to 0(n 3 / log 2 n) for graphs 
with arbitrary edge weights [3]. Obviously, the usage of these methods for 
structures consisting of hundreds of thousands to millions of vertices in many 
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cases is unsuitable or even practically unfeasible due to a large amount of 
computing time. 

Additionally we would like to mention the case where the shortest paths 
between all pairs of a graph's vertices have been calculated while the nu- 
merical characteristics are unknown and are to be found. This happens, for 
example, when initially the computing of numerical characteristics is not 
required or it is necessary to find characteristics not for the whole graph but 
for it's various subgraphs. Moreover, the structures of the subgraphs may 
change over time and, therefore, it is necessary to periodically recompute the 
values of the center, radius and diameter. A trivial solution in this case is 
the examination of all calculated shortest paths, it has a complexity 0(n 2 ). 
The solution of this problem on graphs with large number of vertices or in 
cases where multiple recomputing of graph's characteristics is required may 
also have a long computing time. 

In this paper, two metric characteristics problems are considered. Prob- 
lem 1: is to find the center, radius and diameter of a graph defined only 
by the set of vertices V and the set of edges E and Problem 2: is to find 
the center, radius and diameter of a graph using additional information in 
the form of the calculated shortest paths between all pairs of the graph's 
vertices. 

The problems are considered on connected weighted undirected graphs 
with non-negative real edge weights. Unconnected graphs are not considered 
because, by definition, the radius, center and diameter for them are infinite. 
In cases where characteristics for unconnected graphs are defined as some 
function of characteristics of its connected components (e.g., maximum), 
the problem can be reduced to the considered one by finding the connected 
components with the usage of fast existing methods [9]. Only graphs with 
non-negative edge weights are considered since any graph without negative- 
weight cycles can be associated with the graph having only non-negative 
edge weights with preservation of the shortest paths [8]. 

In this paper, the algorithms for a quick search of metric characteristics 
of graphs for two problems mentioned above are presented. The main ad- 
vantage of the proposed algorithms is that it is not necessary to examine 
all vertices of a graph to determine it's metric characteristics. The pro- 
posed algorithms have shown a significant reduction in computation time 
on weighted graphs of a different nature compared with popular methods in 
solving the problems considered. 

1 Problem definition 

A connected undirected weighted graph G = (V, E, w) is considered, 
here V = (v±, i>2, . . . , v n ) is the set of vertices of the graph, E = (e\, 
. . . , e n ) is the set of edges of the graph, C V x V and w : E — > is a 
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nonnegative real weighted function on edges. It is further assumed that the 
powers of the vertex and edge sets are equal to |V| = n, \E\ = m. Let us 
introduce some notations and definitions. 

Denoted by w(i,j) is the weight of the edge from vertex Vi to vertex 
Vj. The length of the shortest path from vertex vi to vertex Vj is called 
the distance from Vi to Vj and is denoted by rriij. A matrix consisting 
of the shortest paths between all pairs of vertices of a graph is called the 
distance matrix and is denoted by M = (rriij). A graph is called connected, 
if there is a path connecting any two of its vertices, that is the distance 
between any pair of vertices is finite rriij < oo. A weighted graph is called 
undirected, if the weights of the edges between any pair of vertices are equal 
in both directions w(i,j) = w(j,i),Vi,j. Metric characteristics of a graph 
are defined as follows. 

Definition 1. The eccentricity s(vi) of vertex Vi is the maximum distance 
from Vi to any other vertex in the graph, e(vi) = max m^- . 

Definition 2. The radius of a graph r is the minimum eccentricity of any 
vertex in the graph r = mm i= Y^£(vi). Vertex c achieving this minimum is 
called a central vertex (the center), c : e(c) = r. 

Definition 3. The diameter of a graph d is the longest distance between 
any pair of vertices in the graph d = max^ J= x^ rriij ■ Peripheral vertices 
are vertices of a graph such that distances between them are equal to the 
diameter. 

In general, a graph can have several central vertices and several pairs 
of peripheral vertices. Particularly, in a complete graph with identical edge 
weights each vertex will simultaneously be central and peripheral. Our aim 
is to find at least one center and one pair of peripheral vertices, because 
more than often, the knowledge of the radius and diameter is more valuable 
than the information about the vertices on which they are achieved. In this 
paper, two graph metric characteristics problems are to be considered. 

Problem 1. Given a connected undirected weighted graph G = (V,E,w) 
with a nonnegative real weight function on edges w : E — > Find the 

radius r, diameter d and at least one center c and a pair of peripheral vertices 
of the given graph. 

Problem 2. A connected undirected weighted graph G = (V, E, w) with a 
nonnegative real weight function on edges w : E —> K + '° is given and the 
distance matrix M is calculated. Find the radius r, diameter d and at least 
one center c and a pair of peripheral vertices of the given graph. 

Let us introduce two more notions. The single source shortest paths 
problem (SSSP) is a problem in which the shortest paths from a given vertex 
Vi (the source) to all other vertices in the graph vj £ V are to be found. By 
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SSSP(uj) the search of the shortest paths from a vertex Vi to all other vertices 
of a given graph is denoted. The set of pivotal vertices is a set of p vertices of 
a given graph P = (pi,P2, ■ ■ ■ ,P P ) such that the information about distances 
from vertices of this set is used for computing metric characteristics. 

2 Related works 

In a case where the distance matrix of a graph is unknown, the most 
common method of searching for the graph's metric characteristics is to solve 
the all-pairs shortest path problem (APSP). Most of the algorithms solving 
APSP problem aren't universal and they show a good solution speed for 
graphs only with a certain set of properties. In particular, there are algo- 
rithms for sparse graphs [6], for graphs with bounded integer edge weights 
[11] etc. Although the complexity of solving APSP gradually decreases with 
the invention of new algorithms, this approach for the search of metric char- 
acteristics of weighted graphs still remains unpractical for graphs with a 
large number of vertices . 

There are algorithms that find metric characteristics for graphs having 
a special organization and describing peculiar structures. Among these, 
for example, is the algorithm for finding of the diameter for Small World 
Networks graphs [12] or the algorithm for finding of the center and diameter 
for Benzenoid Systems' graphs [4]. To increase the speed of a solution, these 
algorithms use peculiarities of the respective graphs, and therefore the range 
of their effective application is rigidly restricted. 

One of the types of metric characteristics problems on graphs is the 
approximate calculation of the center and diameter of a graph. In these kind 
of problems, distances between graph vertices are searched with some degree 
of error allowing for reduction of computation time significantly. Among 
such, for example, is the Approximate neighborhood function algorithm [2[ 
or [I]. 

Apart from the search of exact and approximate characteristics of a 
graph, the search for the so-called graph effective characteristics is widely 
spread. The effective radius of a vertex v is the 90th-percentile of all dis- 
tances from v. The effective diameter of a graph is the minimum distance 
within which 90% of all vertices are within reach of each other. Certain 
algorithms for such problems are proposed, for example, in |10j . 

A trivial method of solving Problem 2 is to examine the distances be- 
tween all pairs of a graph vertices and determine those satisfying the def- 
initions of the radius and diameter. The solving of this problem can be 
accelerated if computed distances are stored in special data structures (e.g., 
in the sorted list). Otherwise, the construction of such data structures will 
have a time complexity exceeding the time complexity of the trivial method. 
Other published methods of solving Problem 2 are unknown. 
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3 Fast search algorithms 



In short, the main idea of the proposed algorithms is as follows. First, 
a few pivotal vertices are chosen in a certain way and distances from them 
to all other vertices of a graph are computed. Using this information and 
the definitions of graph metric characteristics, the bounds of the radius, the 
diameter and their estimations on current iteration are found. After that, 
vertices having distances out of the found bounds are discarded. If there 
are vertices left, the number of pivotal vertices is increased and the bounds 
and estimations for the radius and diameter are recomputed. The process 
continues until all vertices are discarded. Current estimations of the radius 
and the diameter are the desired characteristics. 

3.1 The center and radius search 

The algorithms for finding the center and the radius for Problems 1 and 
2 are almost identical. The only difference is in the possibility to use com- 
puted distances from a considered vertex. In Problem 2, already computed 
distances, stored in a column of the distance matrix, are simply examined. 
In Problem 1, these distances are computed by the SSSP algorithm. At the 
heart of the search speed-up lies the determination of the lower 77 and upper 
r u bounds of the graph radius with the usage of pivotal vertices P and the 
examination of center pretenders ?«, until the equality 77 = r u doesn't hold 
true. 

The graph radius r has the following property. If the distances from all 
vertices of a graph to an arbitrary vertex Vi are computed and the minimal 
of them rriji = min fe= y^ fc _^ has been found, the graph radius is inferior 
than the lower bound 77 equal to that minimal value r > 77 = rriji. Indeed, 
since rriij is the minimum distance to v i from all other vertices, the eccen- 
tricity of the center c are not lower than this lower bound 77 = mji < m c i < 
e(c) = r. Moreover, 77 is the lower bound of the radius in the case where it 
is defined as the smallest of the longest distances from all vertices of a graph 
to a subgraph of any p vertices. In this case 77 = m i n i=m( max pGP fnip) and 
77 = maxpzprrijp < max pg pm cp < e(c) = r. 

On the other hand, the graph radius r is bounded from above by the 
eccentricity of any vertex r < r u = e(vj). Thus, the graph radius r lies in 
the range 77 < r < r u and 

if n = r u , then r = n = r u (1) 

Thereby, the graph radius r and one of its centers can be found by iterative 
increasing the radius lower bound 77 and decreasing the radius upper bound 
r u until the equality 77 = r u is held. Evidently, the speed of the convergence 
of 77 and r u to the radius depends on the selection of center pretenders q 
and pivotal vertices. 
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For the fast search of the radius, vertices " — center pretenders Cj need 
to be found with as large as possible 

maxm Cj „ = r\ = min (max mj p ) (2) 
peP 1 i=i, n peP 

and with as small as possible eccentricity e(cj), here P " — the set of pivotal 
vertices. It is also desirable that the number of center pretenders Cj examined 
and vertices in P be as small as possible, since it minimizes the number of 
SSSP solved in Problem 1 and the number of the distance matrix elements 
examined in Problem 2. 

Since the central vertex of many graphs is often situated near the geo- 
metrical center (if it could have been found), a good strategy is to choose 
center pretenders having a similar location. One way to find such a selection 
is to choose vertices Cj having the minimum longest distance to peripheral 
vertices of a graph Cj : max^g^ m C! ,j = min i= j-^(maxrf g £) mid), where D is 
the set of peripheral vertices of the graph. This approach requires you to find 
the diameter and peripheral vertices of a graph and therefore is very com- 
putationally complex. However, one can use the following simple method to 
get a good approximation of peripheral vertices. 

The most remote from each other vertices p\ and p2 are sought. 

Pl,P2 ■ m pi P2 = maxm ipi = maxm ip2 (3) 
i=l, n i=l, n 

In order to do that, an arbitrary vertex d\ and a vertex c?2 defined as the 
most remote from d\ di : m,d x d% = max i=TH m dii are chosen. Subsequent 
vertices in this jj sequence,^ are defined similarly 

dj+i : m djd]+1 = maxm^j (4) 

i=l,n 

The process continues until the condition dj—i = dj+i is held, in this case 
the desirable vertices have been found and p\ = dj, p2 = dj + \. 

Vertices p\ and p2 are included in the set of pivotal vertices P = {p\,P2\- 
Center pretenders q are determined by ([2]) . If for the first center pretender 
ci the condition ri = r u is held, c\ is the center of a graph and r = ri = r u . 
If r u > ri, the most remote from c\ vertex p% : m ClP3 = e{c\) is added to 
the set of pivotal vertices P and the search of C2 is performed by ([2]). The 
process continues until the condition ([I]) is held for some pretender. 

By this approach, the search of the graph center and the radius has sev- 
eral advantages. Firstly, the computation (or examination for Problem 2) 
of the shortest paths takes place during the search, and if the condition ([T]) 
is held early enough, there is no need to perform the rest of the procedure 
for most of the graph vertices. Secondly, it is not necessary to compute 
maximums in ([2]) for every new center pretender, one can store computed 
maximums and after adding a new vertex pi to P, compute only maximums 
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from the stored values and the distances from a new vertex p^. The descrip- 
tion of the algorithm for the both problems is represented in fig.l 

The center and radius search algorithm (Rl, R2) 

Input: a graph G = (V, E, w) 
distance matrix M {only for Problem 2). 
n = 0, r u = oc 

1. The search speed-up 

Solve SSSP(dj) for the vertices di from {only for Problem 1). 
Determine pi,p% by ([3|) and include them into P. 

2. Find a center pretender 
Determine the vertex q and r/ by ([2]). 

3. Verification of the pretender 

If the shortest paths aren't computed for a, solve SSSP(q) {only for 
Problem 1). 

Find the eccentricity e(cj) and, if e(cj) < r u , r u = e(cj). 
If r u 7^ ri, add the vertex pj : m CiPj = e{ci) to P and go to the step 2. 
Otherwise, q is one of the graph centers, r = n is the graph radius, the 
algorithm ends. 

Fig 1: The center and radius search algorithm. Additional data and actions 
for Problems 1 and 2 are indicated in italics. 

The time complexity of the algorithm, evidently, depends on the number 
of center pretenders Cj examined. For graphs with different structures this 
number can vary widely, so it is impossible to rigorously evaluate the com- 
plexity for the case in general, nevertheless it is possible to roughly indicate 
the lower and upper bounds. Consider the best and the worst scenarios for 
Problem 1. In the best case, when the first center pretender c\ is one of 
the graph centers, one has to solve the SSSP problem for k vertices from 
the first step of the algorithm, determine c\ by ([2]) and find the eccentricity 
of ci by solving SSSP(ci). The time complexity in this case is determined 
by k ■ dif f{SSSP) + k • 0(n), where dif f{SSSP) is the time complexity 
in solving SSSP. In the theoretically worst case, when the radius is de- 
termined through the examination of all possible pretenders Cj, we obtain 
n ■ dif f {SSSP) + 0{n 2 ), here 0{n 2 ) is the complexity of the eccentricities 
search of all n vertices and determination Cj by ([2]). 

3.2 The diameter search 

The principal difference in the search for the diameter is that, unlike 
for the radius, it is impossible in general to determine a convenient upper 
bound for the diameter. This implies another difference from the radius 
search " — the existence of two different algorithms for Problems 1 and 2. 
Without an upper bound for the diameter the number of pairs of vertices 
examined Vi,Vj " — pretenders on graph peripheral vertices " — can be large 
enough. Hence the examination of each vertex in Problem 1 requires a 
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preliminary solution of SSSP(wfc). Therefore, for Problem 1 it is reasonable 
to perform an auxiliary procedure allowing to reduce the number of pairs 
of vertices Vi,Vj examined and having a smaller computational complexity 
compared to the solving of SSSP. The sorting of vertices pairs vi,Vj by a 
potential distance between them and the examination of the sorted pairs in 
descending order are chosen in the capacity of such procedure. In Problem 
2, the sorting is impractical due to the high complexity of the sorting of 
an array in comparison with the examination of a few extra columns of the 
distance matrix. 

Obviously, the graph diameter d is bounded from below by the distance 
between any two of its vertices d > d\ = rriij,\/i,j = l,n. The upper bound 
of a graph can not be determined without the examination of distances 
greater than the lower bound m^i > d\. That is, the upper bound can be 
computed only after the determination of the graph diameter and thus has 
no practical use. In such circumstances, one way to speed-up the search of 
the diameter is to find a large lower bound di with a minimal computational 
cost and examine only pairs of vertices Vk , vi having distances between them 
potentially exceeding the lower bound m^i > d\. 

The maximal distance between the pivotal vertices is used in the al- 
gorithms as the first approximation of the diameter's lower bound d[. As 
mentioned before, the pivotal vertices are quite a good approximation to 
graph peripheral vertices. The distances from these vertices to the graph 
center c are used in the algorithms for the estimation of potential distances. 
The choosing of c in that role allows us to discard from the examination a 
large number of vertices Vk,vi which are not peripheral niki < dj. 

The diameter search algorithms consist of two parts. At first, one of the 
graph centers is scanned by the algorithms described in the previous section. 
Then the direct search of the graph's diameter, based on the information 
about distances from the center to other vertices of the graph is performed. 
That is, the diameter search algorithms completely solve both considered 
problems and thus may be called the radius and diameter search algorithms. 

3.2.1 Algorithm for Problem 2 

After one of the graph centers c has been found, only the distances from 
vertices Vi for which the following condition is held are examined 



where di is the current lower bound of the diameter determined by 



and P is the set of pivotal vertices generated during the computation of the 
radius. Indeed, if the graph diameter is greater than the current lower bound 



m ic > di/2 



(5) 



di = max ma 
ijeP 



(6) 
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and Vi is one of graph peripheral vertices, for some vertex Vj by the triangle 
inequality for the distance matrix, we have mj c + m c j > rriij = d > d\. 
That is, rriic + m c j > di = 2 • di/2 and it is necessary either mj c > di/2 or 
mcj > di/2. This proves the necessity to consider only vertices Vi satisfying 
([5]), by virtue of a symmetry of the matrix M. The description of the 
algorithm is shown on figure 2. 

The diameter search algorithm for Problem 2 (D2) 

The input: a graph G = (V,E,w), the distance matrix M. 
di = 

1. The search of the graph center 

Find one of the graph centers c by using the algorithm from section 3.1. 
Compute the lower bound of the diameter di by (JBJ). 

2. The search of the diameter 

Examine the distances from vertices Vi satisfying ([5]) may be the diam- 
eter. If rriij > di, d[ = rriij. 
The output: di is the graph diameter. 
Fig 2: The diameter search algorithm for Problem 2. 

3.2.2 Algorithm for Problem 1 

In Problem 1 the examination of all vertices V{ satisfying ([5]) can be quite 
long due to the necessity to solve SSSP(uj) for each such vertex. Therefore, 
to narrow the range of peripheral vertex pretenders one can use the infor- 
mation involving distances from two vertices to the graph center. 

By definition of the shortest path, elements of the distance matrix M 
satisfy the triangle inequality 



Therefore, after having computed distances from an arbitrary vertex v v to 
all other vertices, the distances only between pairs of vertices v^, v\ satisfying 
the following condition have to be examined 



For the vertices not satisfying this inequality, rriki < m kv + m v \ < di, i.e. 
the distances between them are not greater than the current lower bound. 
As the vertex v v the graph center c is used in the algorithm. 

According to (J7J) , vertices vp. , vi with the maximal sum + m v i hav- 
ing the maximal potential distance. Therefore, it is reasonable to examine 
vertices , v\ in descending order of the value m^ v + m v i . If during this 
"top-down" passage for some pair of vertices Vk,vi inequality ((HJ) isn't held, 
the examination of the remaining vertices with smaller potential distances 
between them is not needed. 



mij < rriik + m kj ,\/i, j,k 



(7) 



m kv + m v i > d t 



(8) 
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The sorting of n 2 pairs Vk, vi can be quite long compared to the total time 
of the diameter search, thus one can use the following method which allows 
reduction of the computation time while somewhat degrading the quality 
of the sort. Instead of the sorting of pairs v^,Vi by the sum m^ c + m c i 
one should sort vertices v% by the distance to the center m, c and perform 
verification of condition (jSJ) for vertices pairs in the following order 

"mi i "mi ^m3 j • • • j "mi "m n i "m2 V m3 , !)m 2 ^m4 ■ ■ ■ (9) 

Here m, is the number of the vertex with the ith largest distance to the 
graph center. If some v mi v m . are not satisfying the pair v mi+1 v mi+2 is 
examined. If (|5J) isn't held and j = i + the algorithm is terminated and di 
is the graph diameter. The description of the algorithm is shown in figure 
3. 

The diameter search algorithm for Problem 1 (Dl) 

The input: a graph G = (V, E, w). 
di = 

1. The search of the graph center 

Find one of the graph centers c using the algorithm from section 3.1. 
Compute the lower bound of the diameter d\ by fl5J). 

2. Sorting 

Sort vertices V{ by the distance to the center c in descending order. 

3. The examination of pretenders 
Choose a pair of pretenders from list © 

If pair Vk , vi satisfies (JSj) , then if the shortest paths for Vk , v\ are not 

computed, solve SSSP(-Ufc) and SSSP(u;). 

If m k i > d h di = m kl 

The output: di is the graph diameter. 
Fig 3: The diameter search algorithm for Problem 1. 

As in the case with the search of the radius, the complexity of the algo- 
rithm heavily depends on the graph structure and on how fortunately the 
pivotal vertices were chosen. In general, the complexity of the algorithm 
can be estimated by diff(Ar) + dif f(Sort(n)) + k ■ diff(SSSP) + 0{x). 
Where diff(Ar) is the complexity of the radius and the center search al- 
gorithm, dif f{Sort{n)) is the complexity of the used sorting algorithm, 
k ■ dif f(SSSP) is the complexity of the solving of SSSP for k vertices in 
step 3 of the algorithm and 0(x) is the examination of x pair of vertices 
Vk , v\ to satisfy inequality (JSJ) . 

4 Tests results 

All the tests have been performed on a computer equipped with Intel 
Core 15 530 (2,93 GHz) CPU and 3 GBs of RAM on the 32-bit edition 
of Windows XP. The source code has been written on C++ programming 
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language in Borland C++ Builder 6 IDE. The first set of the data used is 
weighted graphs of the USA road networks from the public source |13j . The 
second set of the data used is complete graphs with random edge weights 
generated by us. For complete graphs results are taken as an average over 
10 runs on different graphs with the same dimension. The parameters of 
the tested graphs are presented in table 1. 



Table 1: The parameters of the tested graphs. 



Name 


Vertices 


Edges 


Average degree 


RN_1 


1001 


2432 


2,43 


RN_2 


2007 


5288 


2,63 


RN_5 


5000 


1,34 • 10 4 


2,69 


RN_7 


7000 


1,8 • 10 4 


2,58 


RN_10 


10 4 


2, 7 • 10 4 


2,69 


RN_15 


1,5 • 10 4 


4, 38 • 10 4 


2,92 


RN_20 


2 ■ 10 4 


5,41 • 10 4 


2,71 


RN_NY 


2,64 • 10 5 


7, 34 • 10 5 


2,78 


RN_BAY 


3,21 • 10 5 


8 ■ 10 5 


2,49 


RN.COL 


4,36 • 10 5 


1,06 • 10 6 


2,43 


RN_FLA 


1,07- 10° 


2, 71 • 10° 


2,53 


RN_NW 


1,21 • 10 s 


2, 84 • 10 6 


2,35 


RN.LKS 


2, 76 • 10 6 


6,89 • 10 6 


2,5 


RN_E 


3,6 • 10 e 


8, 78 • 10 6 


2,44 


RN_W 


6, 26 • 10 6 


1,52 • 10 Y 


2,43 


F_l 


1000 


10 s 


1000 


F_2 


2000 


4- 10 6 


2000 


F_5 


5000 


2,5 • 10' 


5000 


F_7 


7000 


4,9 • 10' 


7000 


F_10 


10 4 


10 8 


10 4 



The proposed algorithms for Problem 1 (denoted by Rl and Dl) have 
been compared to the a) Dijkstra algorithm performed for each vertex of 
graph in the implementation with a binary heap [5] for graphs with a small 
average vertex degree and b) to the Floyd- Warshall algorithm [7] for com- 
plete graphs (denoted by RC1 and DC1). As it is impossible within reason- 
able time to obtain the solution of the problem on high-dimensional graphs 
by using the Dijkstra algorithm, the running time for graphs from RN_BAY 
onward has been approximated using the least squares method. The data 
used for the approximation are subgraphs of different dimensions of the 
tested graphs, a part of them is represented in table Q] (RN_1, RN_2 etc.). 
The running time of the Dijkstra algorithm for all vertices of a graph has 
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been approximated by 7, 1 • 10~ 7 n 2 log(n) — 1,1- 10~ 7 nmlog(n) + 0, 66 for the 
search of the radius and by 7, 2-10~ 7 n 2 log(n) — 1, 1- 10~ 7 nmlog(n)+0, 98 for 
the search of the diameter. The Dijkstra algorithm in the implementation 
with a binary heap has been used for the solving of SSSP in the proposed 
algorithms. Quicksort has been used as the sorting algorithm in Dl. The 
results of the tests for Problem 1 are shown in table [2j 

On all the tested graphs a significant reduction in search time of the 
radius and diameter for Problem 1 is observed when the proposed algorithms 
are used. The percentage of vertices with solved SSSP doesn't exceed 0,9% 
for the radius search and 7% for the diameter search. 

For Problem 2 the tests have been performed on graphs with the number 
of vertices not exceeding 10 4 due to the infeasibility of storing the distances 
information in a computer memory for larger graphs. The proposed algo- 
rithms (R2, D2) have been compared to the simple examination of all ele- 
ments of the distance matrices of graphs for the search of the radius (RC2) 
and to the examination of the upper triangle of distance matrices for the 
search of the diameter (DC2). Due to a small solution time for the problem 
on the tested graphs, each algorithm has been run 1,000 times consecutively. 
The tests results for Problem 2 are shown in table El 

For problem 2 a significant advantage of the proposed algorithms is also 
observed " — the speed-up of the radius search from 56 times onward, the 
speed-up of the diameter search from 5 times onward. 

Conclusion 

The proposed algorithms solve the stated metric characteristics Prob- 
lems 1 and 2 using the distance information only of a small fraction of graph 
vertices. The tests results have shown that, at least on the tested graphs, 
the percentage of the used information decreases with increas of dimensions 
of graphs. This behavior of the algorithms indicates that high-dimensional 
graphs have greater gaining time. 

The main results of the performed tests in terms of numbers are follows. 
For the search of the center and radius for Problem 1 , the number of vertices 
with solved SSSP doesn't exceed 0,9%, the running time speed-up changes 
in the range from 55 up to 1, 28 • 10 6 times faster. For the diameter search 
for Problem 1 these figures are, accordingly, 7% and from 14 up to 1, 22 ■ 10 6 
times faster. For Problem 2 the range of speed-up changes from 56 up to 
1, 14 • 10 3 for the search of the center and radius and from 5 up to 157 times 
faster for the diameter search. 

In further research on the topic of this paper, the test of the proposed 
algorithms on other sorts of graphs and estimation of the sphere of applica- 
tion of the algorithms can be fulfilled. Another direction of research is the 
further optimization of the algorithms through the usage of faster sorting 
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algorithms and the usage of algorithms of solving the SSSP relevant for con- 
sidered graphs. The question of applicability of the approach used for the 
algorithms in search of all graph peripheral vertices and centers also seems 
interesting. 
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Table 2: The tests results for problem 1. Rl SSSP, pieces is the number of the SSSP solving during Rl, Rl SSSP, share is the 
share of vertices with solved SSSP of the total number of the graph vertices. Similarly for Dl. Approximate results denoted 
by*. 



Graph 


RC1, sec 


Rl, sec 


Rl 


Rl SSSP, 


DC1, sec 


Dl, sec 


Dl 


Dl SSSP 








S SSP 7 


share 








S SSP 7 


, share 








pieces 










pieces 




RN_1 


1,33 


0,016 


4 


0,004 


1,31 





4 


0,004 


RN_2 


5,45 


0,031 


7 


0,0035 


5,45 


0,031 


10 


0,005 


RN_5 


39 


0,031 


4 


0,0008 


39 


0,61 


72 


0,0144 


RN_7 


75 


0,063 


6 


0,00086 


75 


5,5 


492 


0,0702 


RN_10 


170 


0,17 


10 


0,001 


170 


0,72 


41 


0,00,1 


RN_15 


373 


0,17 


7 


0,00047 


374 


0,77 


30 


0,002 


RN_20 


761 


0,23 


6 


0,0003 


762 


25 


635 


0,0318 


RN_NY 


1, 54 • 10 5 


3,5 


6 


2,3 • 10- 


-b 


1,54 • 10 5 


5,3 


9 


3,4 • 10" 5 


RN_BAY 


2,49 • 10 5 * 


3,16 


4 


1,3 • 10" 


-5 


2,49 • 10 5 * 


3,31 


4 


1, 25 • 10 -5 


RN.COL 


4, 76 • 10 5 * 


7,09 


6 


1,4-10- 


-5 


4, 77 • 10 5 * 


233 


198 


0,00045 


RN.FLA 


2,99 • 10 y * 


13 


4 


3,7-10" 


-6 


2, 99 • 10 b * 


16 


5 


4,7-10" b 


RN_NW 


4,03 • 10 s * 


14 


4 


3,3 • 10" 


-6 


4, 03 • 10 y * 


15 


4 


3,3 • 10" b 


RN.LKS 


2, 14 • 10 y * 


37 


4 


1,5 • 10" 


-6 


2, 14 • 10 y * 


38 


4 


1,5 • 10" b 


RN_E 


3, 76 • 10 Y * 


49 


4 


1,1 • 10" 


-6 


3, 77 • 10 Y * 


51 


4 


1,1 • 10" b 


RN_W 


1,18 • 10 s * 


92 


4 


6,4-10" 


-7 


1,19 • 10 s * 


97 


4 


6,4- 10~ v 


F_l 


10 


0,18 


9 


0,009 


10 


0,27 


13,9 


0,0139 


F_2 


81 


1 


10,3 


0,0052 


81 


1,75 


17,8 


0,0089 


F_5 


1246 


12 


9,5 


0,0019 


1246 


21 


15,4 


0,0031 


F_7 


3401 


33 


10 


0,0014 


3401 


45 


13,7 


0,002 


F_10 


9853 


71 


9,5 


0,00095 


9853 


123 


16,5 


0,0017 



Table 3: The tests results for Problem 2. 



Graph 


RC2, sec 


R2, sec 


R, speed-up, times 


DC2, sec 


D2, sec 


D, speed-up 


RN_1 


4,21 


0,031 


135 


1,7 


0,047 


36 


RN_2 


156 


0,14 


1114 


7,3 


0,17 


42 


RN_5 


95 


0,16 


593 


45 


1,38 


32 


RN_7 


181 


0,38 


476 


87 


16 


5 


RN_10 


374 


0,89 


420 


177 


2,19 


80 


F_l 


4,07 


0,072 


56 


2,08 


0,107 


19 


F_2 


16 


0,16 


100 


8,16 


0,26 


31 


F_5 


100 


0,43 


232 


50 


0,69 


72 


F_7 


196 


0,63 


311 


98 


0,86 


113 


F_10 


399 


0,81 


492 


200 


1,27 


157 



